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SECTION  I 


INTRODUCTION 


Conventional  techniques  of  simulating  hypervelocity  flight  environ- 
ments in  ground-based  facilities  are  severely  limited  by  their  ability 
to  produce  simultaneously  the  total  enthalpy  and  total  pressure  required 
for  adequate  simulation.  The  testing  mediums  generally  are  not  repre- 
sentative of  true  atmospheric  conditions  because  of  contamination  caused 
by  the  heating  process  and  incomplete  chemical  reccmbination  of  the 
test  gas  dioringjthe.  acceleration  process.  The  multiccmponent  flow  — 
device  described  by  Johnson  and  von  Ohain^  has  the  objective  of 
producing  air  flows  which  simulate  hypervelocity  flight  without  passing 
the  working  gas  through  phases  of  high  static  temperature.  In  this 
device,  a multicanponent  flow  process  is  employed  ii  which  kinetic 
energy  is  transferred  from  a low  molecular  weight  working  gas  to  a much 
higher  molecular  /eight  gas,  such  as  air.  In  both  the  direct  and  in- 
direct energy  exchange  processes,^  either  liquid  air  or  refractory 
particles  are  accelerated  by  the  low  molecular  weight  gas  in  a multi- 
stage process. 

In  addition  to  the  RHE*i  device  of  Ref.  1,  there  is  a critical 
requirement  for  development  of  a combined  erosion-ablation  test  facility 
to  simulate  certain  phases  of  high  speed  vehicle  flight.  In  this 
facility,  solid  particles  are  injected  into  a high  enthalpy  flow;  how- 
ever, the  injection  must  be  accomplished  in  the  supersonic  portions  of 
the  flow  CO  avoid  vaporization  of  the  particles  before  they  reach  the 
test  section.  Accompanying  the  injection  of  peurticles  is  a complicated 
interaction  between  the  main  high  enthalpy  flow  and  the  flow  containing 
the  particles. 

There  is  a paucity  of  data  regarding  focusing  of  solid  and/or 
liquid  particles  in  a two-phase  flow.  The  majority  of  investigations 
on  two-phase  flows  has  been  directed  toward  the  determination  of  the 
effect  of  the  particles  on  the  overall  flow  properties  and  the  deter- 
mination of  particle  trajectories  in  conventional  nozzle  geometries."'"" 
While  the  motivations  for  these  studies  are  much  different  from  those 
of  this  study,  the  techniques  of  analysis  have  direct  application.  The 
theoretical  methods  can  be  used  to  characterize  the  flow  properties, 
and  the  experimental  results  provide  guidelines  for  the  choice  of  par- 
ticle sizes  which  lead  to  optimum  velocity  slip.  When  wave  systems  are 
used  for  focusing,  the  inertial  properties  of  the  particles  and  the 
associated  velocity  and  thenneil  slip  between  them  and  the  carrier  gas 
are  extremely  important  in  determining  the  final  particle  trajectories. 
While  certain  theoretical  analyses  of  gaseous  focusing  systems  are 
available, detailed  experimental  studies  of  focusing  in  two-phase 
flows  are  required  before  successful  operation  of  the  multicomponent 
flow  device  can  be  achieved. 
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A critical  problem  In  the  mtiltlcomponent  flow  device  is  the 
separation  of  the  accelerated  particles  trca  the  low  molecular  weight 
carrier  gas.  In  the  direct  process,  the  particles  must  be  focused  into 
a region  suitable  for  their  vaporization  by  energy  addition.  In  the 
indirect  process,  the  focusing  must  be  acocnplished  in  a manner  which 
leads  to  an  efficient  transfei  of  monentun  and  energy  to  the  tunnel 
working  gas  (air).  Detailed  examinations  of  the  wave  pattern  and  flow 
field  geometries  which  result  in  such  focusing  is  the  subject  of  this 
study. 


Theoretical' and  experimental  studies  have  been  conducted  to  examine 
particle  trajectories  in  various  flow  field  geometries.  Detailed 
analyses  of  shear  drag  coefficients  at  high  Mach  number  and  low  Reynolds 
number  were  conducted,  and  a new  correlation  for  shear  drag  coefficient 
which  Improves  the  accuracy  of  trajectory  calculations  in  high  speed 
flows  was  developed.  Exploratory  experimental  studies  were  conducted 
to  examine  particle  trajectories  in  basic  flow  field  geometries. 
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SECTION  II 


lUO-FHASE  FLOW  THEORETICAL  ANALYSIS 


A.  WO -PHASE  FLOW  CONSERVATION 
EQUATIONS /COMPITTER  PROORAM 

A digittl  oonputtr  progran  haa  baan  davalopad  to  analyca  tha  ataady 
flow  of  a gaa-partlola  nixtura  in  whloh  tha  gaa  affaota  tha  partlolaa, 
but  tha  partialaa  do  not  intaraot  with.tha  gaai  Enargy  axohanga  la 
aaaumad  to  ooour  batwaan  tha  gat  and  tha  partlolaa  by  forced  oonvaotlon 
and  batwaan  tha  partlolaa  and  apace  by  radiation.  Drag  forcaa  ara  tha 
onl/  axtarnal  forcaa  aaaumad  to  act  on  tha  partlolaa.  Thaaa  important 
aaaufflptlona,  in  addition  to  othara  adopted,  ara  llatad  halow. 

1.  Thaoratloal  Tuo-Phaaa  Flow  Aaaumptiona 

a.  In  an  unoouplad  ataady-flow  of  a gaa-partlola  mixtura  vha 
gaa  affaota  tha  partlolaa,  but  tha  partlolaa  do  not  inter- 
act with  tha  gaa. 

b.  No  maaa  or  anargy  la  loat  from  tha  ayataa. 

c.  No  maaa  la  axohangad  batwaan  tha  phaaaa. 

d.  Voluma  occupied  by  tha  partlolaa  la  nagllglbla, 

a.  Thi'mal  (Brownian)  motion  of  partlolaa  la  nagllglbla. 

f.  Partleles  do  not  intaraot. 

g.  uaa  ><nvironmant  la  Inviaold  axoai>t  for  tha  drag  It  axarta 
on  tha  partlolaa. 

h.  Cumpoaitlon  and  apaolflo  haata  of  both  gaa  and  partlolaa 
ara  oonatant. 

I.  Enargy  axohanga  oooura  batwaan  tha  gaa  and  partlola  by 
forced  oonvaotlon,  and  batwaan  tha  partlola  and  apace  by 
radiation. 

J.  Drag  foroaa  ara  tha  only  axtarnal  foroea  acting  on  tha 
partlolaa. 

litlliziti^  tile  above  aaau/njjtluria  tha  conaarvatlpri  aquations  car.  be 
written'-  daaorlDing  the  flow  of  tha  gaa-partlola  ayatam  by  traatlng 
tha  gaa  valooity,  gaa  tamparature,  partlola  velocity,  and  partlola 
tam])arature  aa  dapandant  varlablae  and  dlatanoa  aa  the  Independent 
variable.  Thaaa  aquatlona  have  baan  formulated  In  a coordinate  ayatam 
with  axaa  parallel  and  pari>andloular  to  tha  local  p'lrtlcle  diraotiuii, 
Tha  particle  momentum  equuliona  take  the  form; 


(2) 


In  th«  above  Kuitionf,  the  factor  (fp)  aooountc  for  drvlatlona  of  t'.e 
drag  coeffifllant  from  the  ideal  StoKea  How  value  (24/Ke). 

Similarly,  the  particle  energy  equation  is  given  by 


^ ‘ y,cv»Vla.a)  ‘ -(*P'  - V)]  <3) ■ 


Included  in  the  above  equation  la  a correction  tern  I'g,  alnilar  to 
fp,  to  account  for  deviationa  fren  Stokec  flow,  thereby  including 
rarefaction,  ccopreesibility,  and  inertial  effects.^'’  The  continuity 
equation  and  equation  of  atate  for  a yterfect  gee  complete  the  aet  of 
equationa. 

for  the  purpoabe  of  thia  atudy  it  haa  been  aaaiuned  that  the  gaaeoua 
flow  field  can  be  adequately  deaoribed  by  an  invlacld  two'dinenalonal 
flC''  in  rbich  both  Frandtl' Meyer  expaneion  wavea  and  oblique  ahock  wavea 
can  oe  enbkvlded.  Alao,  the  partlclee  are  aaaiuned  to  be  apherieal  und 
are  dnfined  by  apeclfying  their  molecular  weight,  radiua,  and  apecific 
heat.  There  are  averaging  teohniquea  available  at  MR!  that  are  not 
documentated  in  thia  report  which  allow  mlxturea  of  partlolea  with 
, 'Ifferent  charaoteitatioa  to  be  treated. 

The  viacoaity  tern  which  appeara  in  2q  (l)-(3)  can  be  computed 
frem  the  reaulta  of  Bartt,^  which  la  accurate  for  moat  mixtui'ea  con- 
aiatlng  principally  of  diatomic  gaaea.  However,  for  greater  aocoracy 
the  viacoaity  can  b*  obtained  from  the  reaulta  of  an  analyaia  aimilar 
to  that  of  Adler  and  Anderaon,  * which  la  valid  for  a dilute  gaa  mixture 
with  a frozen  chemical  corniioaitioi.. 

Aa  atated  in  the  initial  aaaumptiona,  the  only  force  acting  to 
accelerate  the  partlolea  ia  the  drag.  The  drag  term  in  Rqe  (1)  and  (2) 
ha:  been  repreaer.ted  by  a ratio  of  the  actual  apherv  dreg  coefficient 
to  that  In  Stokea'  flow,  It  may  be  noted  thet  Stokee'i  flow  regime  ia 
valid  for  Inoompreaalble  flowa  Heynolda  numbara  leaa  than  0.1.  Sinct 
it  ie  poaaiblf  for  the  particle  to  experience  free -molecule,  flip,  and 
continuum  flow  conditiona,  it  la  neoeaaary  to  inoluda  a drat;  oocfflulent 
valid  over  the  entire  rengc  of  flow  oundltlone.  To  thia  end,  the  reeulta 
of  eevercl  drag  coafflolent  oorrelatiune  have  been  uaed  in  thle  atudy 
to  give  em[ilrlcal  flta  of  f.,  ai  a function  of  Mach  ctiu  Heyriolde  nujiibere- 
Theae  are  diacuaaad  In  detail  In  the  next  aeotiun  of  thia  retort, 

A comj/uter  jirogram  to  , mJve  hqa  (l)-(3)  hue  toon 

written  uelng  a I'lmri.h-order  Rmigo-Kut.t  a mimtrricnl  integration  achnmo 
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with  a variable  itep  alze.  A block  diagraa  of  the  method  of  aolution 
ia  given  in  Figure  1.  Caleulatione  begin  at  a apecifled  axial  location, 
Z ■ with  specified  initial  conditions  for  the  gaseous  flow  field, 

i.e.,  Mach  number,  flow  angle,  ratio  of  specific  heats,  total  presstire, 
and  tanperatvre.  The  initial  particle  velocity,  flow  direction,  temper- 
ature, and  concentration  are  specified  as  functions  of  the  radial  coor- 
dinate r.  For  Z > Zaint  selected  particle  trajectories  are  ccmputed 
with  flow  properties  detemlned  from  either  isentropic  expansion  or 
oblique  shock  relationships.  Particle  velocity,  gas  velocity,  particle 
and  gas  flow  directions,  particle  and  gas  teaqperatures,  relative  Mach 
number,  relative  Reynolds  number,  Knudsen  number,  relative  dynamic  pres- 
sure, and  drag  ooeffieient  are  given  at  specified  axial  intervals  and 
corresponding  radial  locations. 

In  the  calculation  of  the  particle  concentration  profiles  the  con- 
tinuity equation  for  the  particle  ia  applied  along  the  previously  com- 
puted particle  trajectories,  which,  In  effect,  are  particle  streamlines. 
Thai  is,  the  continuity  equation  for  the  particles  can  be  written  as 


(cudA)initial  • 


(4) 


where 

, /mass  density  \ /particle\  / . 

^ * \of  particles/ \ velocity / '*’^**' 

Now 

_ /volume  of\  /partlcleX  /number  of  particles\ 

^ “ y particle  / ^density  / \per  unit  volume  / 

i.e. , 

To  find  an  equivalent  form  for  the  udA  term,  consider  two  particle 
streamlines,  1 and  J.  Therefore 


AT  = rj  - r^ 

and 

dA  (Tj  - r^)  • 1 (6) 
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INPUT  DATA 


CALCf  DRAOy  pQ,  Tq,  ^inln>^rnAx 

'r'*' 

[Up,  o,  Tp,  Hp/N^T  - KCNS(r)J 


Fig-jre  1,  Two-Phase  Flow  Computer  Prograri  - Block  Dlagrmn. 


C 


The  velocity  term,  aver-'-ged,  beccmes 


% * 


The  plane  normal  to  the  average  velocity  vector  can  be  expressed  as 


and  the  term  udA.  becomes 


('■j  - '■i 


1 


)L. 


Therefore  the  term  pudA  becomes 


pudA  = 3 nrp^PpNpn  (r ^ ) Icos  j ^ (lO) 

Placing  this  expression  into  Eq  (U)  and  simplifying  result  in 

I®'.  • )i“  ■ 


/initial 


Note  that 
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The  continuity  equation  is  then  satisfied  at  each  specified  axial  loca- 
tion. The  necessary  information  is  available  at  this  location  in  the 
gaseous  flow  field  to  obtain  the  local  particle  concentration  value. 


B.  SPHERE  DRAG  COEFFlCIiNT  CORREIATIONS 

An  important  ingredient  in  the  two-phase  flow  theoretical  analyses 
is  the  correlation  employed  to  obtain  the  drag  coefficient  of  the  solid 
particles.  Various  drag  coefficient  correlations  were  examined  and 
were  compared  with  the  results  of  experL-nental  sphere  drag  coefficient 
measurements  obtained  in  a ballistic  range. These  data  provide 
sphere  drag  coefficients  over  a range  of  Mach  numbers  from  approximately 
0.1  to  7.0  with  sphere  Reynolds  numbers  ranging  from  5 x 10®  to  8 x 10^, 
This  set  of  data  includes  over  500  data  points  taken  in  the  same  facilitv 
by  the  same  investigators.  These  data  are  for  T^/T^  = 1 and  therefore 
comparisons  of  sphere  drag  coefficients  can  be  made  without  considering 
wall  temperature  effects. 

The  sphere  drag  coefficient  correlation  of  Csurlson  and  Hoglund^^ 


2U  r 

" Re  L 


ii. 


0.15 

~TI 


+ (M/ReK3.B2 


expC(O.U27/M~^-^^)  - 

+1.28  exp[-1.25  Rey 


(13) 


and  that  of  Crowe 


Cp  (Cp^^-'')exp[ -3.07(7  )"'(N'/Re)g(Re)]  + exp[ -Re/2M]  + 2 

where 

logiog(Re)  = 1.25(1  + tanh(0.77  logioRe  - 1.92)1  (lU) 


and 


h(M) 


2.3  tanh(l.l7  logiQM) 


were  first  examined.  Comparisons  of  the  drag  coefficients  of  Carlson 
and  Hoglund  and  Crowe  are  shown  in  terms  of  percent  deviation  from 
experimental  data^^  in  Figure  2.  Percent  deviation  is  defined  as 
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Figure  2,  CompELrison  of  Existing  Sphere  Drag  Correlations  with  Elxperimental  Data 


% Deviation  = 


Ct^ 

^theor 

exp 


X 100 


(15) 


As  can  be  seen,  the  differences  between  the  theoretical  and  experimental 
sphere  drag  coefflci^its  of  Carlson  and  Hoglund  are  extensive,  ranging 
to  ± 80^,  while  for  the  correlation  of  Crowe  the  deviations  range  f^on 
+U0^t  to  -30^>  It  is  interesting  to  note  that  the  Carlson  and  Hoglund 
predictions  are  consister<-ly  high  for  experljnental  sphere  drag  coeffi« 
eients  greater  than  ai>proximate)y  1.6.  Conversely,  Crowe  predicts  both 
high  and  low-theoretical  sphere  drag  coefficients  over-  the  entire  rcuige 
of  experimental  sphere  drag  coefficients.  A significant  difference  not 
apparent  in  this  comparison  is  that  the  drag  coefficient  correlation  of 
Crowe  allows  for  various  particle  to  gas  temperature  ratios  while  that 
of  Carlson  and  Hoglund  includes  no  wall  temperature  effects. 

Another  sphere  drag  coefficient  correlation  that  was  examined  is 
that  of  Cuddthy,  Beckwith,  and  Schroeder.^' 

For  M < 0.5  1 


Cp  = 


51.1M 

Ree 


1 . o.256m(c„^  . 


where 

exp^-0.028  Reg' 


For  M > 0.5 


Cd  = (^D,FM  ' ^D,c)  [-ARes"^] 


where 

0.520  < 
2.00  < 
O.OUU7  < 
0.410  < 


Cd,C  < 0*981  I 
0d,FM  S 7 -So  1 
A < 0.315  I 
n < 0.745  I 


(Table  Inputs  as 


FNS  of  Mach  No.) 


This  correlation^^  comparea  favorably  with  the  ballistic  range  data  of 
Bailey  and  Hyatt, as  shown  in  Figure  3,  with  maxiJBum  deviaticxis  rang- 
ing from  approximately  +305t  to  -bo^.  It  does  not,  however,  include  the 
influences  of  unequal  particle  and  gas  tonperatures. 

A new  sphere  drag  correlation  formula  has  been  developed  from  the 
ballistic  range  data  by  modifying  the  correlation  of  Cuddihy,  Beckwith, 
and  Schroeder  to  include  varying  particle  and  gas  temperatures. 

For  K < 0.*)  1 


Cd  = 


51. IM 
Re- 


1+0 


.256m 


51. IM^ 

Re„  > 


where 


CDq  = ^ + 0.U4  + 1.6  exp[-0.028  Res®**""'] 


For  M > 0.5 

Cd  = h3.(M)[l  + exp(-ARes’^)]  - hr>(M)exp (-ARCj.")  + <££!d(T) 


(17) 


where 


0.455  < hi(M)  < 1.045  \ 

2.020  < h.(M)  < 7.800 

0.045  < A < 0.315 


5.410  < n < 0.740  | 


} (Table  Inputs  as  FTJS  of  Mach  No.) 


Tw 

1 


^d(t)  - “7^  + 2.220  - 0.597(^p]  exp[-ReVM]|  (M  > l) 


Comparisons  of  Eq  (l6)  and  (17)  show  that  for  M < 0.5  the  correlation 
is  identical  to  that  of  Cuddjhy,  Beckwith  and  Schroeder.  However,  for 
M > 0.5  the  rearrangement  of  the  C^  expression  and  the  changes  to  the 
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table  inputs  of  hi.(M}  and  h,o(M)  result  in  a more  accurate  fit  of  the 
experijMntal  data  (sae  Table  I for  new  table  inputs).  A comparison  of 
the  percent  deviation  of  the  results  of  this  present  correlation  with 
the  ballistic  range  data  is  shown  in  Figure  3>  It  can  be  seen  that  the 
present  correlation  greatly  improves  the  accuracy  of  the  sphere  drag 
coefficients  t resulting  in  deviations  of  to  -20^  over  an  extremely 
wide  range  of  slip  Mach  number  and  slip  Reyn^ds  number.  These  com- 
parisons have  been  made  for  Ty/Tj,  =1.  As  shown  in  the  present  corre- 
lation (Eq  (17)],  a ^|)(T)  tern  has  been  added  to  include  the  effects 
of  variable  wall  temperatures.  The  results  of  Bailey  and  Hiatt have 
again  been  u&el  to  verify  the  accuracy  of  the  ^j)(T)  term,  as  shown  in 


Figures  U throxigh  7 
the  freestream  slip 
normal  shock  (Reo). 
ccmparlson : 


It  may  be  noted  that  the  correlation  is  based  on 
Mach  number  (M)  and  a slip  Reynolds  number  behind 
The  following  cases  have  been  used  for  means  of 


a 


M = 2.0 

M = U.O 
M = 6.0 

M = 10. U 


Reo  = 6,10,20,50,100 

Reo  = 2,7,10,20,50,200,700 
Re-  = 7,10,20,50,200,1000 

Re-  = 2,5,10,20,50,500,5000 


(Figure  b) 
(Figure  5) 
(Figure  6) 
(Figure  7) 


As  shown  in  Figures  4 through  7,  the  <y!Jj)(T)  term  prediction  is 
acceptable  for  a range  of  Mach  numbers  from  2 to  10.4  and  Reynolds 
numbers  frcm  2 to  5000.  At  the  present  time,  the  <5£:])(T)  term  is  re- 
stricted to  M > 1 since  no  comparisons  have  been  made  for  the  subsonic 
slip  Mach  numbers. 


C.  INFLUENCE  OF  SPHERE  DRAG  COEFFICIENT  CORRELATION 
ON  TWO-ffiASE  FLOW  CHARACTERISTICS 

The  Influence  of  the  sphere  drag  coefficient  correlations  on 
particle  trajectories,  velocities,  and  temperatures  in  simple  flow 
situations  was  investigated.  The  flow  situations  consisted  of  uniform 
flow,  where  the  particles  are  injected  at  an  angle  from  the  horizontal 
with  a specific  initial  velocity,  a simple  Prandtl-Meyer  expansion  fan, 
and  an  oblique  shock  wave.  The  cases  investigated  were  for  particle 
radii  of  0.5,  5>  25,  and  50  um  and  for  particle  velocities  of  100,  1000, 
and  2000  ft/s.  In  this  segment  of  the  study,  titanium  dioxide  (Ti02) 
particles  were  used  and  the  fluid  was  taken  as  air  (molecular  weight  = 

29). 


1.  Uniform  Flow 

The  trajectories  for  a 1.0  um  diameter  particle  injected  at  an 
initial  velocity  of  100  ft/s  into  a Mach  5 air  flow  at  an  angle  of  45* 
are  shown  in  Figure  8.  Note  that  there  are  large  differences  in  the 
particle  trajectories,  depending  on  the  specific  drag  coefficient  corre- 
lation employed.  Hence,  accurate  particle  drag  coefficients  are  re- 
quired to  obtain  reliable  particle  trajectories.  Corresponding 
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TABLE  I 


TABLE  INPUTS  AS  FUNCTIONS  OF  MACH  NIMBEK 
FOR  FRESOTT  SPHERE  DRAG  CORRELATION 


Mach  No. 

hi(M) 

hs(M) 

A 

n 

0.5 

0.1+80 

7.800 

0.315 

0.410 

,0.6 

0.U55- 

6.500 

0.240 

0.468 

0.7 

O.U97 

5.570 

0.182 

0.500 

0.8 

0.535 

4.920 

0.141 

0.545 

0.9 

0.610 

4.450 

0.110 

0.590 

1.0 

0.820 

4.100 

0.090 

0.620 

1.1 

0.900 

3.850 

0.070 

0.645 

1.2 

0.990 

3.600 

0.065 

0.670 

1.3 

0.995 

3.420 

0.060 

0.680 

l.U 

0.995 

3.230 

0.055 

0.690 

1.5 

0.999 

3.110 

0.052 

0.700 

1.6 

1.000 

2.980 

0.049 

0.710 

1.7 

1.020 

2.890 

0.048 

0.713 

1.8 

1.C05 

2.800 

0.047 

0.715 

1.9 

0.99C 

2.740 

0.047 

0.713 

2.0 

1.045 

2.680 

0.046 

0.720 

2.2 

1.035 

2.580 

0.046 

0.723 

2.U 

1.040 

2.480 

0.046 

0.725 

2.6 

1.035 

2.420 

0.046 

0.725 

2.8 

1.035 

2.360 

0.046 

0.725 

3.0 

1.035 

2.320 

0.045 

0.727 

3.2 

1.030 

2.280 

0.045 

0.730 

u.o 

1.025 

2.170 

0.045 

0.730 

5.0 

1.030 

2.100 

0.045 

0.735 

6.0 

1.040 

2.050 

0.045 

0.735 

8.0 

1.030 

2.C20 

0.045 

0.740 

IK 


"he  Compeu-ison  of  Tmperature  Dependent  Sphere  Drag  Data  with  the  Present 
Correlation  (M^  = 2,0). 


iture  Dependent  Splieare  Dra^  Data  with  the  Present 


variRtions  of  the  particle  velocity  and  particle  temperature  obtained 
with  the  v6u:iou8  drag  coefficients  are  shown  l,i  Figures  9 and  10.  It 
can  be  seen  that  there  is  little  influence  of  the  specific  drag  coef- 
ficient correlation  on  the  particle  velocity  and  temperature. 

Since  the  gas  velocity  in  this  case  is  232U  ft/s,  an  initial 
particle  velocity  of  2000  ft/s  was  chosen  to  investigate  the  effect  of 
this  variable.  The  results  are  shown  in  Figure  11,  where  again  it  is 
noted  that  large  differences  in  the  particle  trajectories  depending  on 
the  specific  drag  coefficient  correlation  are  found.  Also,  as  shown  in 
Figures  12  and  13  the  corresponding  variations  of  peo-ticle  velocity  and 
psurticle-temperature  show  little  influence  of- the  drag  coefficient 
correlation . 

The  particle  size  in  the  previous  cases  was  relatively  small. 
Therefore,  the  case  where  a 100  urn  diameter  particle  was  injected  in 
air  at  an  initial  velocity  of  100  ft/s  at  an  angle  of  U5®  was  investi- 
gated. The  particle  trajectories  for  the  sphere  drag  coefficient  cor- 
relations for  this  case  are  shown  in  Figure  14.  It  is  interesting  to 
note  that  for  this  particle  size,  there  is  no  significant  difference  in 
trajectories,  as  experienced  with  the  1 mUi  particle  diameter,  for  the 
four  different  sphere  drag  correlations.  However,  as  shown  in  Figures 
15  and  16,  there  are  significant  differences  in  particle  velocity  and 
temperatures.  Hence,  accurate  sphere  drag  coefficients  are  required 
for  all  particle  sizes  to  correctly  predict  the  trajectory,  velocity, 
and  temperature  of  the  particles. 

Other  cases  were  investigated  for  particle  injection  into  a 
uniform  flow  and  are  contained  in  Appendix  A.  These  cases  arc  for 
peirticle  radii  of  0.5,  5,  25,  and  50  ..m  and  for  initial  particle 
velocities  of  100,  1000,  and  2000  ft/s.  The  results  of  these  catlcula- 
tlons  are  similar  to  those  discussed  above. 

2.  Prandtl-Meyer  Exi)ansion  Fan 

A series  of  numerical  experiments  were  conducted  to  examine 
particle  trajectories  through  simple  expansion  waves.  Typical  results 
for  a 1 urn  diameter  particle  with  an  initial  velocity  of  100  ft/s 
traveling  through  a I randtl-Meyer  expansion  fan  with  a total  turning 
angle  of  30"  are  shown  in  Figure  17.  Results  are  given  for  all  four 
drag  coefficient  correlations  examined,  and  it  is  clecr  that  the  drag 
coefficient  correlation  has  a significant  effect  on  the  particle  tra- 
jectories through  the  expansion  fan.  It  is  also  interesting  to  note 
that  the  corresponding  particle  velocity  and  particle  temperature  vari- 
ations through  the  expansion  fan  shown  in  Figures  18  and  19  show  little 
difference  in  the  sphere  drag  coefficient  correlation  used. 

As  in  the  unlfomi  flow  investigation,  this  f;ase  was  also 
examined  for  an  initial  velocity  of  2000  ft/s.  The  resulting  particle 
trajectories  s!'0wti  in  Figure  20  display  even  more  significant  differ- 
ences. Notably  the  results  with  the  present  sphere  drag  coefficient 
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Figure  10.  Particle  Tenperature  History  in  Air  for  Uniform  Flow  Case  for 

Different  Drag  Coefficient  Correlations  (r_  =0,5  um,  u_  = 100  ft/s) 


Carlson  & Uoglund 

Beckwith,  Cuddlhy  t Schroeder 


Trei^iectory  of  Particle  in  Air  for  Uniform  Flow  Case  for  Different  Drag 
Coefficient  Correlations  (r  = 0.5  nm,  u = 2000  ft/s),' 


cle  Velocity  History  in  Air  for  Uniform  Flow  Case  for  Different 
Coefficient  Correlations  (r„  = 0.5  nJn»  u~,  = 2000  ft/s). 


100  ft/i 


Carlson  4 HoRlund 
Beckwith,  Cnddlhy  i Schroeder 


Figure  l6.  Particle  Temperature  History  in  Air  for  Uniform  Flow  Case  for 

Different  Drag  Coefficient  Correlations  (rp  = 0,5  M®»  = 100  ft/s) 
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correlation  are  similar  to  those  using  the  Cuddlhy,  Beclwith,  and 
Schroeder^^  correlation.  The  corresponding  particle  velocity  and  par> 
tide  temperature  variations  are  shown  in  Figures  21  and  22.  These  a 

results  indicate  a significant  difference  in  the  particle  velocities 
3 to  U feet  dovmstream  from  the  turning  point;  however,  there  is  little 
difference  in  the  particle  temperattires . 

As  a means  of  comparison,  the  case  for  a 100  pm  diameter  par- 
tide  with  an  initial  velocity  of  100  ft/s  is  shown  in  Figure  23.  Here 
again,  the  trajectories  for  the  larger  particle  show  little  effect  due 
to  changes  in  the  sphere  drag  coefficient  as  opposed  to  those  for  the 
1 Mffl  diameter  particles.  The  corresponding  particle  velocities  shown 
in  Figure  2U,  and  the  particle  temperatures  given  in  Figure  23  show  ; 

greater  effects. of  changes  in. the  drag  correlation..  _ _ . . j 

Additional  cases  were  also  investigated  for  the  Prandtl-Meyer 
expansion  fan  and  are  conttined  in  Appendix  B.  These  cases  are  for  j 

particle  radii  of  0.5,  5,  and  50  um  and  for  initial  particle  velocities  j 

of  100,  1000,  and  2000  ft/s  turning  an  angle  of  30* . ] 

3.  Two-Dimensional  Wedge/Obllque  Shock 

Numerical  experiments  were  cond\  - ted  tc  determine  the  trajec- 
tories, velocities,  and  temperatures  of  the  particles  when  passing 
through  a two-dimensional  oblique  shock  wave.  A 30*  wedge  in  a Mach  5 
air  flow  which  generated  an  oblique  shock  with  an  angle  of  approximately 
42*  was  chosen.  The  trajectory  of  a particle  with  a 10  urn  diameter  and 
an  initial  velocity  of  100  ft/s  was  computed  for  the  four  sphere  drag 
coefficient  correlations,  with  the  results  shown  in  Figure  26.  It  can 
be  seen  that  there  is  no  significant  difference  in  the  particle  tra- 
jectories. This  type  of  behavior  is  also  exhibited  in  the  particle 
velocity  and  particle  temperature  variation  as  shown  in  Figures  27  and 
28.  Therefore,  in  the  compression  case  there  is  little  variation  in 
results  as  the  sphere  drag  coefficient  is  changed.  This  is  substantiated 
by  the  additional  cases  examined;  these  results  are  Included  in  Appendix  ^ 

C for  particle  radii  of  5 and  50  um  and  for  initial  particle  velocities 
of  100,  1000,  and  2000  ft/s.  1 


D.  COMPOSITE  SOLWION  TO  DESCRIBE  FHE  PARTICLE 
FLOW  ABOUT  A TWC-DIMhWSIONAL  WEDGE 


The  previous  section  described  the  particle  behavior  in  simple 
flow  situations,  i.e.,  xmlform  flow,  Prandtl-Meyer  expansion  fan,  and 
an  oblique  shock.  The  flow  field  generated  by  a wedge  when  placed  in 
a supersonic  stream  would  contain  all  three  of  these  flow  situations. 
Therefore,  a 10*  half-angle  two-dimensional  wedge  was  chosen  and  placed 
into  a Mach  2 air  flow,  and  the  computer  program  was  used  to  determine 
the  two-pliaae  Inviscid  solution.  This  calculation  was  performed  in  two 
steps.  First,  given  the  free  stream  conditions  and  the  wedge  half- 
angle, one  obtains  the  solution  for  a two-dimensional  infinite  wedge  in 
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a supersonic  stream.  This  solution  yields  the  flow  variables  and  the 
particle  trajectories  at  all  points  behind  the  shock. 

To  ccoplete  the  solution  for  the  finite  wedge,  the  expansion  of 
the  flow  about  the  shoulder  of  the  wedge  defined  by  the  point  (ry,z^) 
must  be  computed.  To  obtain  this  solution,  the  flow  field  variables 
must  be  specified  as  functions  of  r at  the  initial  axial  location, 

These  initial  conditions  are  detemined  from  the  results  of  the  two- 
dimensicxial  infinite  wedge  solution.  It  should  be  noted  that  the 
Prandtl-Meyer  expansion  cannot  be  used  in  its  general  form  for  the  flow 
over  the  shoulder  of  the  wedge  since  theory  predicts  an  expansion  to  an 
infinite  Mach  number  and  zero  pressure.  Physically,  there  is  a free 
shear  layer  formed  at  the  sho\ilder  of  the  wedge  which  must  be  accounted  - 
for  in  the  calciilations  to  obtain  the  appropriate  turning  angle  of  the 
flow.  In  the  present  computer  program,  the  flow  conditions  behind  the 
wedge  are  specified  with  a base  pressure  correlation^”  for  flows  over 
two-dimensional  bodies.  From  this  correlation  the  base  pressure  behind 
the  wedge  can  be  calculated  for  a given  free  stream  Mach  number,  and 
from  this  base  pressure  the  proper  turning  for  the  flow  can  be  calculated. 

The  resiilts  of  the  application  of  this  procedure  to  a two-dimensional 
10°  half-angle  wedge  for  a 2 pm  diameter  titanium  dioxide  particle  are 
shown  in  Figure  29.  The  following  free  stream  conditions  were  used: 

P-^  = 7920  PSF  Up^  = 100  ft/s 

T-  = 540*R  Tp  = 540°R 

-,=14  ^ 

Z % = 0° 

M^  = 2.0  ^i 

arid  results  in  the  particle  trajectories  shown  in  Figure  29.  The 
velocities  and  temperature  variation  for  two  particle  trajectories  are 
shown  in  Figures  30  and  3I.  In  addition,  the  mass  flux  density  profiles 
at  three  axial  locations  in  the  flow  field  are  shown  in  Figure  32  and 
tend  to  emphasize  the  focusing  effect  when  a bo(iy'  of  this  type  is  placed 
into  a supersonic  two-phase  flow.  Hence,  it  is  possible  to  predict  the 
entire  two-phase  flow  field  about  a body  of  arbitrary  symmetry  placed 
into  a supersonic  stream  since  the  flow  field  can  be  constructed 
from  the  uniform  flow,  Frandtl-Meyer  expansion  fan,  and  oblique  shock 
flow  cases. 


E.  EFFECT  OF  DIFFERBUT  GASES  ON 
WO -PHASE  FLOW  CHARACTERISTICS 

Investigations  were  conducted  to  determine  if  anj'  significant 
difference  results  for  particles  injected  into  a gas  different  from  air. 
A uniform  flow  was  examined  for  the  following  gas  conditions: 
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The  ;artlclc6  were  aieumed  to  be  titanium  dioxide  with  a diameter  of 
10  tun,  injected  at  an  angle  of  1*5*  with  an  initial  velocity  of  100  fv/s 
and  an  initial  temperature  of  5^0^R.  The  particle  trajectories  shown  in 
Figure  33  correspond  to  particle  injection  into  helium,  air,  and  argon. 

As  can  be  seen,  the  radial  coordinate  is  directly  proportional  to  the 
molecular  weight  and  there  is  a direct  viscosity .effect.  The  gas  vis* 
oosity  was  assmed  proportional  to  the  molecular  v?elght  of  the  gas.^ 

Since  the  Mach  number  was  specified,  the  gas  velocities  for  the 
three  cases  investigated  were  different.  The  ratio  of  particle  velocity 
to  gas  velocity  is  shown  in  Figure  3^.  Again,  there  is  a definite  re- 
lationship between  the  magnitude  of  Up/ug  and  the  molecular  welgnt. 

ThtLj  results  indicate  that  Up/ug  aj/proaches  unity  quicker  for  the 
larger  molecular  weight  gas.  This  is  sl^iif leant  when  considering  the 
momentum  of  the  particle  available  for  tremsfer.  In  contrast  to  these 
results,  the  ratio  of  particle  temperature  to  gas  temperature  was  used 
to  evaluate  the  effect  of  different  gojes  on  the  energy  exchange  process. 
Ttie  results  are  shown  In  Figure  35.  There  is  little  difference  between 
results  for  argon  and  air*,  however,  Tj,/Tg  is  larger  for  helium  than  it 
is  for  both  argon  and  uir. 

F.  EFFfcCT  OF  DIFFTIdiJiT  I’AKTlCLEH  ON 
W0-I'IIA;;E  FUjW  CHARACTEHIOTIC;; 

In  this  phase  of  the  study,  the  flow  situation  described  in  the 
previous  section  was  udu]ted  fur  nir  and  the  tyjie  of  partlcl''  was  varied : 
i.e.,  Tin  , and  AI  u , [/'iiticlee  were  used.  Tlie  characteristics  of 

these  ]iarticles  ure; 


lartlcle 

Molecular 

V/elght 

Donoi  ^..v 

ft  ') 

6])ecific  Heat  Capacity 
(ft  /f  -“R) 

TKt  . 

79.90 

8.6 

75.64 

W 

40.32 

10. 0 

9680 

Al  0., 

101.90 

7.7 

8100 

'I  ho  renult J til'  p'lrtlnio  trajoctorloo  oro  shown  in  Kl^pjrc  3''  '-ind  indicate 
thero  j ;;  n nlrnl fj cut.'.  dllTuroiici.'  oviui  1 honi'h  thoro  la  u largo  differ- 
once  iii  tho  tnoloculur  wcl/'.li'  o,  'ihio  io  uloo  true  for  jiartlcle  vclocl- 
tJ()S  uiiiJ  t,uj;i)  oriturott , '.n  nhown  li.  Fiipjroo  37  '^Uid  3^'.  i ho  cunclusion 
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SECliOK  III 

TWO-mSE  FLOW  EJCPERIMmAL  ST'JDIES 


An  experimental  study  to  examine  wave  patterns  which  lead  to 
focusing  of  refractory  particles  was  undertaken  in  this  segment  of  the 
study.  A special  wind  tunnel  facility  was  assembled,  and  solid  par- 
ticles of  known  diameter  and  concentration  were  added  to  the  nozzle  flow 
as  a first  order  simulation  of  the  jet  interaction  region  in  a multi- 
con5)onent  flow  device.  Bodies  generating  various  wave  systems  were 
inserted  into  the  two-phase  flow  field,  and  laser  scattering  techniques 
were  employed  for- the  determination  of  local  particle. concentration. 


A . TEST  FACILITIES 

The  experimental  studies  were  conducted  in  a continuous  flow  wind 
tunnel  system  at  the  Aeronautical  and  Astrcnautical  Research  Laboratory. 

The  wind  tunnel  employs  an  air  storage  system  of  1500  ft*  at  2650  psig 
and  a two-stage  ejector  system  located  in  the  diffuser.  The  run  times 
for  the  facility  are  limited  only  by  the  gas  storage  system.  A sche- 
matic drawing  of  the  facility  is  shown  in  Figure  39.  It  consists  of  a 
conical  nozzle  (exit  diameter  of  1.75  inches)  with  an  exit  air  Mach 
number  of  5 exhausting  to  an  open  jet  test  cabin.  The  measured  pitot 
pressure  profile  for  this  facility  taken  at  the  nozzle  exit  plaiie  is 
shown  in  Figure  Uo. 

Particles  were  injected  in  the  stagnation  r>;iion  of  the  facility 
with  a fluidized  bed  system  shown  schematically  in  Figure  4l,  A helium 
supply  was  used  to  pressurize  the  particle  container  and  thus  fluidize 
the  particles  within  the  porous  inner  container.  The  particles  were 
then  injectel  into  the  tunnel  by  opening  the  on-off  valve  contained  in 
the  feed  tube.  The  flow  of  particles  was  further  regulated  with  a 
0,090-iach  orifice  in  the  feed  t be. 

The  particles  were  titanium  dioxide,  manufactured  by  Dow  Chemical 
Company.  A particle  sample  was  analyzed  with  a Coulter-Coujiter  Model  T 
'Jc.  102"^  to  obtain  a:,  accurate  meas’-urement  of  the  particle  size  distri- 
butic.  . The  results  are  shewn  in  Figiore  42  in  terms  of  a particle  size 
histogram,  i.e.,  percent  differential  volume  versus  particle  diameter. 

Percent  differential  volume  is  defined  as  the  percent  of  the  total 
volume  occupied  by  the  particles  in  a specif' c particle  diameter  range; 
e.g.,  approximately  2Un-,  of  the  volume  is  occupied  by  particles  ranging 
between  0.8  and  1.0  .im  in  diameter.  In  addition  to  the  TiO-,  the  par- 
ticle sample  contained  approximately  1'';  by  weight  of  CAVASIL  whose 
diaiaeter  was  on  the  order  of  0,1  .jro.  This  material  when  added  to  the 
TiO  has  a liquidizing  effect,  thereby  increasing  the  flowing  qualities 
of  the  powder.  The  particle  histogram  (Fig,  42)  does  not  reflect  this 
addition.  i 


55 


j 

9 


F 


Figure  39.  Particle  Concentration  by  Wave  Systems  - Facility  Schematic 


pitot  probe 


Measured  Pitot  Pr-jsure  Profile  for  Two-Phase  Flow  Tunnel 
(No  Particles). 


B.  MDDELB 


Varioui  wtdgt  modcla  w«r«  inacrttd  Into  th«  two’phait  flow  to 
dctormlno  th«  offoot  of  tht  gacnotry  and  aMOolatad  wav#  ayattn  on  tha 
partlola  oonoantratlon  profllaa.  Wadgaa,  5 and  19*,  apanning  tha  antlra 
nozzla  axit  plana  with  baaa  dlmanalona  of  0.088  aj'id  0,262  inch,  raapac* 
tlvaly,  wart  uaad.  Tha  modtla  wart  aat  at  0,  2.8,  7.8 i and  10*  anglaa 
of  attack. 


C.  018TRUMWTAT10N 

Vgrtcua  partiela  ooMt&tratim  naMuraaaati  taohaiquaa  hava  baan 
d^lopad,  prlaarlly  for  uac.in  rookat  notora.  Tbaia  Ineluda  apao*  - 
troioopio  aad  raoovary  tank  mathoda,^**  maa  aaapla  prebaa  uaad  with 
maaa  apaetronatara^  and  ohronotoi  .-aphaf  ’ and  p^tograpblo  obaarva- 
tion  taohnlquaa*''  Howavar,  thaaa  mathoda  althar  parturb  tha  flow 
fiald  baoauaa  of  tha  praaanoa  of  a aanpllag  proba  or  giva  intagratad 
raaulta  aoroaa  tha  flw  fiald.  It  la  alao  unlikely  that  tha  gaa 
aaii'pling  uathoda  will  yiald  raliabla  maaauramanta  of  partlola  oon> 
ownbration, 

Kaoantly,  laaar  aoattarlng  tichnlquaa  hava  baan  davalopad  for  tha 
mMoauramant  of  j4i''tlolt  oonoantratlona.  In  theie  taohnlquoa,  a laaar 
baaxn  ia  pro.laotad  aoroaa  tha  flow  fiald,  and  acattariiig  of  tho  photona 
by  aolid  partiol'-a  raiult.  Ti^-i^al  apjdloatlona  of  tha  taohnlqua  can 
ba  found  in  the  litaratura.  Tha  typ^  cf  laaar  aoattarlng  which 

raaulta  daponda  upon  tha  ratio  of  tha  partlola  aiza  to  wavalangth  of 
tha  laaar  light,  V/hen  the  pajtiola  aJza  oomparad  to  the  wavelength  la 
aikull,  Haylaigh  aoattarlng  r«aulta.  'iha  work  of  Daum  and  Karrall'  " for 
tha  dataotlon  of  air  oondanaatiori  in  hyparaonio  wind  tunnal  la  typloel 
of  tha  application  of  Raylalgh  acattarlng,  Ho!'/tv#r,  whan  aolid  par* 
tiolaa  ara  praaant  in  tha  flow  field,  aoattarlng  raaulta  with  a oroaa- 
aaotlon  oonaidarably  largar  than  that  for  tha  Raylalgh  aoattaringi 
tharafora,  tha  atvara  raqulramanta  j>lacad  on  iha  aoattarlng  dataotlon 
ara  ullaviatad  when  solid  |'')rtlalaa  are  roa]>oitalbla  for  tha  Inaor 
■oattaring, 

Laaar  ecattarlnK  tauhnlqua^t  ware  cmployod  in  the  ]>reaant  atudy  f^r 
tha  naaauramant  of  portiola  concoritrutiona . Tha  Inatrumantotlon 
rniitio  for  thia  aysiem  la  ahowr.  in  Plgura  •*'2,  A 3 ^ lialiuin«naon  laser 
waa  uaad  and  tha  aoattojad  radiation  wna  niaaaura'l  vrlth  a ]>hotomultlpliar . 
Thia  ayetam  glvaa  accurate  maaauramania  of  iooal  partial#  oonoantratlona 
at  aalaotad  point!  within  tha  two«phase  flow  field.  Olnoa  the  laser 
beoin  ia  relatively  unaffaotad  by  tha  flow  field,  good  ejjatiai  roaulutlon 
in  tha  maaeuraroanta  la  obtainad. 

Particle  oonoantratlona  ware  ‘napj.ad  by  aendlnc  tha  haom  from  tha 
laaar  througti  tha  teat  cabin  and  meaaurlng  tha  acattared  light  Intenalty. 
;jcuna  wore  made  along  the  laaor  Onm  and  IriteriBlty  voraua  Ulatarico 
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profiles  were  obtained.  The  output  from  the  spectrometer  was  amplified 
and  processed  by  an  analog  computer  and  recorded  on  an  X-Y  plotter. 


D.  RESULTS  AND  DISCUSSION 

A study  was  first  conducted  to  determine  the  characteristics  of 
the  two-phase  flow  with  the  present  injector  scheme.  The  spectrometer 
was  fixed  at  the  geometrical  centerline  of  the  facility,  and  the  par- 
ticles were  injected  into  the  flow.  The  recorded  intensity  magnitude 
as  a function  of  time  is  shovm  in  Figure  44  and  indicates  that  there  is 
an  approximate  three-second  segment  in  which  the  flow  of  particles  is 
essentially  constant.  With  the  present  physical  arrangement,  it  appears 
that  the  0.090-inch  orifice  is  effective  only  during  this  time  span. 
However,  since  the  filter  used  in  this  study  contains  only  a small 
quantity  of  particles  which  are  injected  through  a fluidized  bed. 

Figure  44  may  indicate  the  result  of  a decreasing  number  of  particles  . 
available  for  injection.  This  conclusion  is  verified  by  the  data  of 
Figure  4$,  where  the  percent  decrease  in  intensity  magnitude  was 
measured  at  a fixed  location  as  the  valve  was  opened  and  closed  a number 
of  times  without  refilling  the  particle  cylinder.  The  results  show  that 
the  fluidized  bed  must  be  refilled  for  each  measurement  to  obtain  con- 
sistent data. 

1.  Basic  Particle  Concentration  Profile 

A typical  nondimensional  intensity  profile  is  shown  in  Figure 

46.  The  particle  concentration  profile  is  strongly  peaked  on  the  tunnel 
centerline  when  no  wave  systems  are  present  to  deflect  the  particles. 

The  corresponding  photographs  taken  in  the  facility  are  shown  in  Figures 

47,  48,  and  49.  . 

2.  Particle  Concentration  Profiles 

Behind  Wedges 

Particle  concentration  profiles  downstream  of  various  aero- 
dynamic shapes  have  been  investigated  with  the  laser  scattering  tech- 
nique. Photographs  of  the  scattered  laser  light  downstream  of  a 5 and 
15“  wedge  at  angles  of  attack  of  0,  2.5,  7.5,  and  10°  are  shown  in 
Figures  50  through  55.  The  variation  in  intensity  of  the  scattered 
light  is  clearly  evident  and  reflects  the  variations  of  local  particle 
concentration  downstream  of  the  wedge.  These  results  indicate  the 
focusing  effect  that  can  be  obtained  by  placing  an  aerodynamic  shape 
and  its  associated  wave  system  into  a two-phase  flow. 
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igure  51.  I’rofilfe  of  Titanium  Dioxide  ; articles  Due  to 
Scattered  L.ight  - 5*  Total  ;\ngle  Wedge  at  0" 
Angle  of  Attack. 


SECTION  IV 


SUMMARY  AND  CONCLUSIONS 


The  experimental  and  theoretical  results  of  this  study  have  pro- 
vided new  information  regarding  the  behavior  of  solid  particles  in  high 
speed  flow  systems.  A new  drag  coefficient  correlation  which  allows 
greater  accuracy  in  the  prediction  of  particle  trajectories  over  a wide 
range  of  slip  Mach  number  and  slip  Reynolds  nvunber  was  formtilated.  A 
variety  of  basic  two-phase  flow  situations  were  investigated,  i.e., 
uniform  flow,  Prandtl-Meyer  expansion,  and  oblic^ue  shock,  by  varying 
the  initial  particle  velocity  and  particle  radius.  In  addition,  the 
effects  of  different  gases  and  different  particles  on  the  two-phase 
flow  characteristics  were  investigated  and  significant  differences  were 
noted.  The  exploratory  experimental  studies  conducted  in  this  investi- 
gation substantiate  that  a test  facility  suitable  for  the  study  of  par- 
ticle trajectories  has  been  developed  and  that  accurate  particle 
concentration  profiles  can  be  obtained  with  the  laser  scattering 
technique. 
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UNIFORM  FLOW 
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This  appendix  contains  the  following  cases: 


a 

n 

(um) 

Up 

(ft/s) 

(r  ^ z) 

(Up  ^ z) 

(Tp  vs  z) 

0.5 

100 

45 

0.5 

1000 

X 

45 

0.5 

2000 

45 

5 

100 

X 

X 

X 

45 

5 

1000 

X 

45 

5 

2000 

X 

45 

25 

100 

X 

X 

X 

45 

25 

1000 

X 

45 

25 

2000 

X 

45 

50 

100 

45 

50 

1000 

X 

45 

50 

2000 

X 

Note:  Thfcse  numerical  calculations  were  carried  out  for  the 


following  conditions: 

M',’  = Air  (29) 

P:,  = 7876.8  psf 

"'‘Pi  = ^5° 

To  = 540'’R 

?g  = 0” 

M,  = 5.0 

Tp  ^ 540"R 

Titanium  Dioxide 

•ox 


Crirl;-.0!i  & i'Dijlund 


Ison  t Hoglund 

Ixlch,  Cuddlhy  t Schroeder 


§ 


t 


U O tl 

3 Si  &t 


O G O I 


SW*Ki 


6.0 


Carlson  & Hoglund 
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PRANDTL-^ER  EXPANSION  FAN 


Pnceding  page  blank 


95 


f 


Carlson  & Hoglund 
Beckwith,  Cuddlhy  1 Schroeder 


(FT) 


Oarlsar.  & Hoc'und 


Carlson  k Hoglund 
Beckwith,  Cuddihy  t Schroeder 


o 

tJ 

c 

8 8 

O CL, 


o n o 


« 

a 

o 


1000  ft/t 


Carlson  & Hoglund 
Eeckwlth,  Cuddlhy  & Schroeder 
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APPE3JDIX  C 


TWO-DBIENSIONAL  WEDGE/OBLiQUE  SHOCK 
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This 

appendix  contains  the  foU.owing 

cases : 

a 

n 

(un) 

Up 

(ft/s) 

(r  ^ 

(Up  ^ z] 

1 (Tp  ^ z) 

30 

5 

100 

30 

5 

1000 

X 

30.  , 

_5  . 

2000 

X 

X - 

X - - 

30 

50 

100 

X 

X 

- 

30 

50 

1000 

X 

30 

50 

2000 

X 

Note: 

These  mimerical  calculations  were  carried  out 

for  the 

following  conditions: 

IV 

- Air  (29) 

= 

7876.8  psf 

^i 

= 0* 

T-  = 

540*R 

= 30‘ 

= 

5.0 

"Pi 

= 540* H 

Titanium  Dioxide 
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(FT) 


arisen  4 Hoglund 

ecicwlth,  Cuddlhy  4 Schroedej 


Carls 


(FT) 


(FT) 


LIST  OF  SYMBOLS 

Defined  by  Eqe.  (16)  and  (I7)i  erea 
Specific  heat 

Drag  coefficient  baaed  on  croaa-aectlonal  area 
Specific  heat  at  constant  pressure 

Dlaneter--  - . _ _ . 

Concentration 

Defined  by  Eqs.  (l6)  and  (17) 

Distance  In  the  relative  flow  direction 

Reynolds  number  based  on  slip  condition  behind  a normal  shock, 

Re^  c 

Radius 
Temperature 
Velocity 
Axial  distance 

Particle  flow  angle  measured  from  the  horizontal 

Angle  from  the  horizontal  of  the  relative  velocity  vector 

Ratio  of  specific  heats 

Particle  emiasivity 

Angle  measured  from  the  horizontal 

Viscosity 

Density 

Stefan-Boltanann  constant 
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SubacrlptB 


c Continuum 

FM  Free  molecule 

g Oas 

Inc  Incompressible 

p Particle 

rel  Relative  between  particle  and  gas 

T Total 

w Wall 

2 After  normal  shock 
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